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contents necessarily reflect the views and policies of the 
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INTRODUCTION 


The continental slope of the Gulf of Mexico is poorly known from an 
ecological standpoint, but is a habitat into which oil and gis exploration 
and development activity is proceeding at an accelerating rate. In 
response to its mandate to obtain environmental data on the impacts of 
petroleum development on the ovter continental shelf (OCS) and provide the 
relevant information to decision makers, Minerals Management Service (MMS) 
contracted in late 1983 a four-year research program to investigate the 
Slope habitat across the Gulf of Mexico waters north of 25°N and from 200 
to 2600 m in depth. This study is being conducted by LGL Ecological 
Research Associates, Inc. (LGL) in cooperation with Texas A&M University 
(TAMU). The first year of the program has been completed and is reported 


herein. 


PROGRAM OBJECTIVES 


The overall objectives of the program are: 


(1) To determine the abundance, structure, and distribution 
of animal communities in the deep sea in the Gulf of 
Mexico. 

(2) To determine the hydrographic structure of the water 
column and bottom conditions at selected sites within the 
Study area. 

(3) To determine and compare sedimentary characteristics at 
selected sites within the study area. 

(4) To relate differences in tiological communities to 
hydrographic, sedimentary, and geographic variables. 

(5) To assess seasonal changes in deep sea biological 
communities in terms of abundance, structure, animal 
size, and reproductive state. 

(6) To measure present levels of hydrocarbon contamination in 


the deep-sea sediments and selected animals prior to, and 


in anticipation of, petroleum resource development beyond 
the shelf-slope break. 

(7) To compare the biological and non-biological 
characteristics of the deep Gulf of Mexico with that of 
other temperate and subtropical deep-sea regions. 

(8) To assemble together and synthesize appropriate published 
and unpublished data with the results of this study, 
summarizing on a seasonal and spatial basis all 
biological, habitat, and other environmental observations 
and parameters. 

(9) To conduct an effective quality assurance and quality 
control program which insures that all data acquired are 
accurate and repeatable within standards normally 
required for each type of observation, measurement, or 
determination. 

(10) To critically review, interpret, and analyze all 
observations and data acquired to redefine as necessary 
the research program in such a way as to avoid or 
minimize redundancy and to optimize the efficiency of all 
field, labortory, and data management operations for 
future deep-sea studies sponsored by MMS in the Gulf of 
Mexico. 

(11) To assess the need for and determine the type of studies 
to be conducted in future program efforts. 


The data for meeting these objectives have been gathered on five 
cruises of which two were conducted during this first year of the progran. 


STUDY AREA AND METHODS 


During Year One, three transects of five stations each were sampled, 
with one transect located in each of the Western, Central, and Eastern 
Gulf of Mexico Lease Planning Areas (Fig. 1). Average sampling depths 
along each transect were 348, 657, 839, 1341, and 2530 m which 
corresponded to the approximate bathymetric centers of previously 
hypothesized faunal zones: Shelf/Slope Transition (150-450 m); 
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Figure 1. Location of transects and stations, within Western (W), Centi.' (C) and Eastern (E) Gulf of Mexico lease 
planning areas. 


Archibenthal, Horizon A (475-750 m); Archibenthal, Horizon B (775-950 m); 
the Upper Abyssal (975-2250); and the Mesoabyssal, Horizon C (2275-2700 
m). 

The types and numbers of samples taken at each station included 800 
bottom photographs; a hydrocast and supportive w " column sampling, 
either six (Central Transect) or three (East and West Transects) box 
cores; and a trawl sample. Box core samples were split for biological 
(macroinfauna, meiofauna) and physical (sediment grain size and chemistry) 
analyses. In November 1983, only the Central Transect was sampled whereas 
in April 198%, all three transects were scmpled. All samples planned to 
be taken were obtained. 

The laboratory analyses are proceeding on schedule using standard 
protocols for each type of physical, chemical, and biological analysis. 
The results which were available are reported below and in the main report 


volume. 


ENVIRONMENTAL FEATURES 


The water column over each transect was characterized by the presence 
of distinctive water masses (Fig. 2) that showed little seasonal or 
geographic variation. From top to bottom, these water masses included a 
shallow mixed layer of Gulf Water (usually present from the surface to 250 
m); Tropical Atlantic Central Water (“300 to 500 m); Antarctic 
Intermediate Water (“500 to 1000 m); and Gulf Deep Water, a mixture of 
North Atlantic Deep and Caribbean Mid-water. The Shelf/Slope Transition 
Faunal Zone roughly corresponds in depth to the Tropical Atlantic Central 
Water mass, the Archibenthal Faunal Zone to the Antarctic Intermediate 
Water mass, and both the Upper Abyssal ani Mesoabyssal faunal zones ocour 
in Gulf Deep Water. Gulf Deep Water is distinctly colder than the water 
masses above, and temperature is undoubtably one of the important factors 
controlling depth distributions of organisms. 

Bottom sediments at stations in the Central Transect in November 1983 
(Fig. 3) were predominantly clay-sized at Stations C1, C2, and C3, with 
deeper stations containing higher proportions of either silt ([5) or siit- 
and sand-sized particles (C4). The same stations on the Central Transect 
during Cruise II (April 1984) typically had a higher proportion of silt- 
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Figure 2. Water masses along the Central Transect during Cruise I. A similar distribution 


was observed for Cruise II. 
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Figure 3. Sediment grain sizes at Cruise I sampling stations. 
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sized particles than had been observed in November of the previous fall 
(Fig. 3 and Fig. 4%). Sediment levels of organic carbon (Fig. 5) and 
calcium carbonate (Fig. 6) were also higher in samples taken on Cruise II 
than in samples takenonCruise I. Results of the sediment hydrocarbon 
levels also suggested that an influx of terrigenous material (bulk organic 
matter and plant biowaxes) to the bottom occurred between the two 
samplings. 

On a geographic basis, bottom sediments at stations on the Eastern 
Transect contained considerably more sand and silt than sediments on the 
other transects, even though all stations were predominantly clay. 
Calcium carbonate levels were highest in sediments from the Eastern 
Transect, and higher in Western Transect samples than in Central Transect 
samples. The pattern of organic carbon levels in the sediments indicated 
levels to have been highest on the Central Transect, and then generally 
higher for sediments from the Western Transect than for sediments from the 
Eastern Transect. An exception was noted for the deepest station; i.e., 
organic carbon levels at Station E5 were higher than levels at Station WS. 
Organic carbon levels exhibited a trend of decrease with depth. At the 
Central Transect, sediment organic carbon in Nevewber 1963 was 
characteristic of carbon provided by marine phytoplankton, based upon 
carbon isotopic analyses. 

With one exception, results from carbon isotopic analyses for benthic 
organisms not collected in the vicinity of gas seeps in April 1983 
Suggested that the slope biota derive their energy from sinking 
photosynthetic carbon (marine phytoplankton). The exception (a crab, 
Geryon guinguedens) had a carbon isotopic value suggesting a food source 
other than marine phytoplankton alone. Animals collected from around 
seeps had carbon isotope levels suggesting chemosynthesis, as opposed to 
photosynthesis, provided the primary production at the base of the seep- 
community food chain. 

Sediments at all three transects had a mixture of thermogenic, 
terrigenous, and planktonic hydrocarbons. The two samplings at the 
Central Transect suggested an influx of low UCM terrigenous material 
occurred between Cruises I and II. This terrigenous material consisted 
primarily of bulk organic matter and plant biowaxes. The material being 
transported to this area appeared to be compositionally constant with 
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time. The biowaxes were accompanied by a low molecular weight UCM and by 
n-C,< to n-Cyq compounds. The higher molecular weight UCM appeared to 
accumulate in place and was much more highly degraded than the terrigenous 
material. Piston coring in the Gulf of Mexico intraslope has demonstrated 
that the Central Transect is in an area of active natural oil seepage. 
Piston cores sampled at these sites generally showed an increase in 
hydrocarbons with depth. This suggests that the source of the high 
molecular weight UCM in the sediments is upward migration, though 
transport of anthropogenic hydrocarbons to the sediment by water column 
particulates cannot be ruled out. 

The influence of riverborne material in the sediments decress< from 
the Central to the West to the East Transect. The reduced hydrocarbon 
levels in the East Transect were primarily due to smaller terrigenous and 
thermogenic inputs. Planktonic and algal inputs were difficult to discern 
in the West and Central Transects, but were readily apparent in the East 
Transect as shown by the numerous alkenes detected. This may be due to 
more rapid sedimentation rates at the Central and West Transects and/or 
the large input of riverine material causing rapid dilution of oceanic 
detritus. Elevated microbial activity in the sediments and/or in the 
water column may also assist in removing the more labile marine debris. 

In general, hydrocarbons were only present in low concentrations in 
the sediments, especially at the East Transect. Aliphatic hydrocarbon 
levels ranged from “10 to 50 ppm. Aliphatic hydrocarbon levels recorded 
in the literature range from 1 to 3000 ppm. The low concentrations 
generally occur in very sandy areas, whereas the high concentrations occur 
in polluted, shallow waters. In areas of pervasive seepage on the Gulf of 
Mexico slope, aliphatic hydrocarbons have been measured in excess of 
100,000 ppm. 

With one exception, all organisms surveyed for hydrocarbon 
contamination appeared to be pristine. The exception was a pooled sample 
of shrimp (Nematocarcinus rotundus, five individuals) from Station E3. A 
complete suite of alkanes and the unresolved complex mixture present in 
this sample strongly suggested petroleum contamination. However, bottom 
tars were also collected in this trawl. The shrimp may have become 
contaminated in the trawl, but one would expect that contamination during 
Sampling would have been confined to the exterior hard parts. 
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BIOLOGICAL FEATURES 


Biological studies include investigations of the meiofauna (organisms 
passing through a 0.30 mm screen but retained on a 62 micron screen), the 
macroinfauna (organisms retained on a 0.30 mm screen) and the megafauna 
which were sampled by trawling. The meiofaunal collections from the 
Central Transect for Cruises I and II indicated a substantial increase in 
density occurred in April 1984 as compared to levels observed in November 
1983, particularly at the shallowest station (Fig. 7). There was also a 
marked change in the relative abundance of major taxa, namely ‘he 
increased relative abundance of Foraminifera. The density data for the 
Central Transect suggested a trend of decreasing abundance of meiofauna 
with depth. 

Comparisons of the three geographic regions in terms of meiofauna 
density and composition showed the Central Transect to have had the 
highest levels as well as a higher proportion of Foraminifera than was 
present in collections from the other two transects (Fig. 8). Whereas the 
collections from the Western Transect exhibited a moderate decline in 
abundance with depth, no such trend was observed for samples from the 
Eastern Transect. 

The meiofaunal collections are yielding good numbers of a newly 
described phylum, the Loricifera, and the poorly-known kinorhynchs. Both 
groups contain at least genera and species new to science. It is now 
known that the Loricifera are not restricted to shallow, sandy substrates, 
given that collections have been taken as deep as 2530 m on clay bottoms. 

Density levels of macroinfauna from the slope taken during this study 
are markedly higher (range was from 2435 to 8628 organisms/m2) than levels 
previously reported from the Gulf of Mexico slope and abyss (25 to 1095 
organisms/m@), Our samples were screened with a seive having 0.30 mm mesh 
whereas the previous study used a 0.42 mm mesh. The typical macroinfauna 
we are seeing from the samples are minute, making weighing impractical 
without destroying the samples. The size feature probably also accounts, 
in large part, for the disparity between our and previous measurements of 
macroinfaunal density from deep-Gulf habitats. 

The seasonal cata for the macroinfauna from the Central Transect also 


suggested an increase in density in April 1984 as compared to November 
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Figure 7. Comparison of meiofauna densities (no./10 cm-) between Cruises I 
(November 1983) and II (April 1984) on the Central Transect. Note 


the large increase of forams in Cruise II diagrams. Comparisons can 
be made on both axes. 
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Figure 8 . Comparison of meiofauna densities (no./cm?) obtained during 
Cruise II (April 1984). Comparisons can be made on both axes. 
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1984, but the increase was not nearly as pronounced as the change observed 
for meiofauna (Fig. 9). On the Central Transect, density of macroinfauna 
did not exhibit a pronounced decrease with depth, but there was an obvious 
decline in abundance at 2530 m (Station ©5). Polychaetes were the 
numerical dominants at all depths on the Central Transect excert at 2530 mf 
where nematodes of a macroinfaunal size equalled or exceeded the relative 
numbers of polychaetes. 

Macroinfaunal densities on the Eastern and Central Transects were 
higher than the density of organisms found on the Western Transect (Fig. 
10). The Western Transect also differed in that macroinfaunal density 
exhibited a decline with depth. On both Central and Eastern Transects, 
density levels were rather consistent from 348 to 1341 m, dropping sharply 
at 2530 a. 

The macroinfaunal groups which have been sorted to the species level 
are exceedingly diverse and contain many new species and genera. The 
taxonomist for the tanaidacean collections has indicated that once all the 
specimens from Cruise II are described, it will increase the number of 
known species, world-wide, by 20%. 

Of the megafauna collected to date, the decapod crustaceans, 
echinoderms and demersal fish collections have been identified to the 
species level. There were 78 species of decapods (led in variety by the 
anomurans, including the galatheids), 33 species of echinoderms (not 
including the brittle stars) and 94 species of fish representing 42 
families. Bathymetric distributional patterns of the megafauna collected 
to date agree very closely with previous work, providing credence to 
historical faunal zonation and assemblage characterization schemes based 
upon the megafauna. 

The benthic photography aspect of the program to date has been mainly 
devoted to the development of quantitative analytical procedures which 
have now been finalized. Preliminary results of photographic analysis for 
Stations Wi, C1, and E1 from Cruise II indicate Station C1 was 
characterized by a greater density of both biota and Lebensspuren than the 


other two stations. 
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Figure 10. 


Comparison of macroinfaunal densities (no. /m?) obtained during 
Cruise II (April 1984). 
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CONCLUSIONS 


Few conclusions can, or should, be made at this early point in the 
program. However, it would appear that there arc marked regional 
differences in the slope environment and biota, as weil as seasonal 
changes. The latter may prove to be related, in large part, to the 
influence of river discharge. 

From the standpoint of hydrocarbon contamination, the slope 
environment and biota have thus far appeared pristine, or nearly so. 
Natura! seeps are prevalent in the vicinity of the Central Transect and 
may, in fact, provide an additional source of energy to newly discovered 
chemosyntnetic communities in this region. Data from Cruise III suggest 
this to be the case. The results of that cruise will be described in 


subsequent reports. 
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ABSTRACT 


This study involved the modification and subsequent application of a 
generic plume model developed by the Walden Division of Abcor Inc. under con- 
tract to the US Army Corps of Engineers. The mathematical basis of the model 
is a unidirectional, steady, Reynolds type diffusion equation. To close the 
governing equation, coefficients of eddy diffusion are introduced and a hypo- 
thesis of similar mass and momentum diffusion invoked to define them. Also 
made is an assumption that the velocity may be replaced by its mean over the 
depth. This substitution facilitates separation of the original equation into 
two partial differential equations. Ome equation models the effect of lateral 
diffusion and is solved analytically. The other models the interaction of 
settling, longitudinal transport, and vertical diffusion and is solved numer- 
ically. The effect of this separation of variables is a substantial reduction 
ir the computational labor. 

As a part of the study a new computer voroeram was wrirten tn nerfarm the 
numerical computations. The reasons for this were 


1) The original model only allowed for initial, upstream distributions 
which were uniform over the depth. The new model allows for inflow 
of the material through a "window" which begins at some elevation 
at or above the bottom and ends at or below the free-surface. 


el required that the distribution of material at the 
y be specified as concentration. The revised model 
ncentration or mass flux to be used. 


2) The original mod 
upstream boundar 
allows either co 


3) The original model used (apparently) a general library subroutine 
for solving systems of numerical equations. It was desired to have 
a version of the program which could be used reasonably on much 
slower microcomputers. Therefore, a semi-implicit, tridiagonal algo- 
rithm was introduced to make the computations highly efficient. 


4) Output tables not generated by the original model were desired and 
code appropriate for generating these was added. Principal among 
these is a table which gives the rate of deposition. 


A second program which sets up the data needed by the numerical compu- 
tation program was also written. This program asks the modeler a series of 
questions and stores the responses in a data file which is subsequently used 
by the numerical plume model. The combination of the two programs form an 


interactive, easy to use, tool for studying mud plumes. 
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PART I: INTRODUCTION 


l. The study which is described in this report involved the modification 
of a generic plume model to compute the dispersion of material discharged in 
off-shore drilling operations. The original model was developed by 8B. A. 
Wechsler and D. aA. Cogley* under contract to the Corps of Engineers as a par* 
of the Dredged Material Research Program. The model assumes that the velocity 
field is time-invariant and constant in magnitude and direction. The depth of 
the fluid is also constant. The planer region occupied by the fluid is as- 
sumed to be a half-plane. This latter assumption is consistent with the envi- 
ronment of off-shore drilling rigs which are usually far from shoreline 
boundaries. 
tie vojectives of tne study were: 


a. Make the model more flexible so that cases where the initial 


distribution of sediment varied with respect to depth could be 


b. Incorporate a nonconstant definition of eddv diffusion into the 


ncorporate an algorithm into the model to comput 
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sediment flux at various downstream distances. 


Incorporate an algorithm into the model to compute the rate of 
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deposition. 

e. Adapt the model for execution on small micro- or mini-computers. 
These are in essence the tasks which were accomplished. Relative to item (e), 
an examination of the model was made in terms of its numerical efficiency. 
It was determined that the efficiency of the computations could be much in- 
proved by taking advantage of the fact that the system of equations to be 
solved was tridiagonal and stationary so that the equation needed to be solved 
only once instead of for each progressive step downstream from the injection 


point. | 


* Wechsler, B. A. and Cogley, D. R., “A Laboratory Study of the Turbidity 
Generation Potential of Sediments to be Dredged,’ WES Technical Report 
D-77*14, November 1977. 
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Emperical relations are required to close the analytical description 


of a turbulent-plume. Wechsler and Cogley use an algebraic closure which is 


not unreasonable, although, conversely it is not necessarily correct. Pursuit 


of this point is beyond the scope of the current investigation. 
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PART II: DESCRIPTION OF THE MODEL 


Basic Equations 


3. The basic equation to be solved is the steady state three-dimensional 


transport-diffusion equation. Assuming eddy diffusion, this equation can be 


written as: 


3 3 ‘ ; 3c 3 3c > f~ 3¢ 
— (uc) + — IJwf(w)dw- — (E —)- — (E =) - = [e. =): 0 (1) 
3x Vv x K ox vy y oy 12 Z 32 
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downstream vertical eddy 
advection sedimentation diffusion 


where 
x = downstream coordinate 


vertical coordinate 
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z = lateral coordinate 

u ® current velocity at any point 

c * sediment concentration 

w ® settling velocity 

f(w) = settling-velocity frequency distribution 

EE ye. = eddy diffusivities in x, vy, and z directions 

4. This equation is based on the assumptions that the flow is steady, 
uniform, and fully turbulent, and that eddy diffusion can be characterized by 
Fick's law with eddy diffusion coefficients. Equatior | is unclosed and sev- 
eral assumptions were made by Cogley and Wechsler in order to simplify its 


solution. These are: 


a. Eddy diffusion in the downstrean. direction is negligible com- 
pared to the other two diffusive transport terms. 

b. The flow is fully turbulent, therefore, the local velocity u 
in Equation 1 can be replaced by the mean velocity U. 

c. The mass and momentum diffusion coefficients in the vertical 


direction are essentially the same. 


d. The shear stress is linearly distributed over the depth, and 


the velocity distribution is given by Von Karman's 


“universal” distribution. 


e. The lateral eddy diffusion coefficient is essentially constant 
over the depth. 
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5 Based on the above assumptions, algebraic expressions for E and 


E. can be written as: 


E = 0.4 U y(l - y/h) sqrt(f/8) (2) 


and 
E = 0.2 hU sqrt(f/8) (3) 


where f is the Darcy friction factor and U is the depth averaged current 


velocity. For each fraction of the sediment load for which the settling veloc- 


ity is constant, Equation 1 may be rewritten as: 


aC C3 ne Cg (4) 


v= saee as £ os == : 
ox sv ay Vv av 12 Zz 92 


rocess considered 


anal 
rv 
as 
® 

wo 


which together with Equations 2 and 3 is the model o 


nere. 


Boundary Conditions 


6. Physically, the boundary condition for the free surface must specify 
that there is no fiux across the free surface. This condition is described 


mathematically by: 


Cogiey and Wechsler assume that all of the sediment reaching the bottom is 
deposited and that there is no reentrainment. The corresponding boundarv 


condition is: 


Separation of Variables 


7. It may be assumed that the solution for Equation 4 is of the form: 


Ci{x,y,z) * Co (x,y) C4 (x52) (7) 


UJ 


On substituting Equation 7 into Equation 4 the following two equations are 


obtained: 


aC aC, ; aC 
U—-+vw c— ii <= (8) 
x dy av \y dy 
and 
oo 2 ( “2) (9) 
x Iz \z 32 


Or, stated another way, if C, and C, are solutions of Equations 8 and 9, 
respectively, their product is a solution of Equation 4. A solution of Equa- 
tion 9, which is applicable to situations where there are no lateral con- 


:traints on the flow is: 


| -1/2 h 
47xE (2-v)“u 
) = = a a 
C(x, 2) ; i exp XE, dv (10) 


Equation 8 is a mathematical description of a two-dimensional diffusion process 


aa 


in which there is no lateral diffusion. It must be solved numerically. 

8. The plume's suspended sediment concentration i.e. the solution given 
by Equation 7, is symmetric about the x-axis and the maximum concentration occurs 
at the centerline. Equation 9 mcdels the distribution of the sediment in the 


1 -eral direction. The quantity C, (x,0) is then in effect a description of 


d 
tne centerline dilution which results from lateral diffusion. 

9. In concluding this section it seems pertinent to observe that Equa- 
tion 4 is homogenous which means that if C(x,v,z) is some solution then so is 
aC(x,y,z) where a is any constant. The consequence of this ,omogenuity is 
that C has no particular dimensions associated with it. The convention which 
has been adopted is to specify distances in meters, velocities in meters/sec, 
and concentration in grams per cubic meter (equal to milligrams per liter). 


It follows that mass fluxes have the dimensions of grams/sec and that rates of 


deposition have the dimensions of grams/sec per square-meter. 
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Numerical Approximations 


10. A numerical approximation of Equation 8 


Cartesian grid. The approximation emploved int 


in this section. Equation 8 may be rewritten as: 
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li. Using a Taylor series expansion it may be shown that to an order of 
dy squared: 


2g Sw (CE (jl) 48, (§))/(2¢dy*dy) *C_ (4, 541) 
3 j y 


+ (CE (Cj) +E (j-1))/ (24dy*dy)*C_ (4, j-1) 


- CE (j+1)+24E (3)+E (j-1))/(2*dy*dy)*C (i,j) (15) 
This relationship is used as a replacement for che left hand side of Equa- 
ction 11 in the numerical approximation. 

2. The result of substituting Equations 14 and 15 into Equation 11 is 
a system of linear equations for each i. Examining Equation 15, it can be 
seen that the system of equations is tridiagonal and furthermore that the sys- 
tem is independent of i, i.e., x. These factors influence the tise re- 
quired to accomplish the computation since the system of equations need be 
solved only once. Alsc, while in general a computation of order n cubed is 
required to solve n linear equations, the computation required to solve a 
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tridiagcnal system is only of order jn. 


PART III: APPLICATION OF THE MODEL 


Model Input 


13. The input quantities which must be specified are the depth h , the 
mean current velocity U , the Darcy friction factor f , and a description of 
the sediment and the manner ‘n which it enters the flow. The sediment is re- 
garded as consisting of a finite sum of parts or fractions, and the fall velo- 
city for each fraction must be specified. The sediment enters the water column 
through an imaginary or virtual window (see Figure 1). The modeler, having 
specified the depth, h , specifies the number of stations in the y or verti- 
cal direction, Nys. Station Number 1 is at the bottom and the station Nys is at 


the free surface. The distance between stations in the vertical direction is 


dh = h/(Nys - 1) (16) 


The bottom and top of the window are specified in terms of station numbers. 


For example, if 


h = 50 meters 


and 


then 


and if the bottom of the window is specified as vertical station 10 and the 
top of the window is specified as station 15 then 

elevation o; window bottom = 10*2.5 = 25 meters 
and 

elevation of top of window = 15*2.5 = 37.5 meters 
Sediment inflow into the mcdel may be specified in terms of either a mass con- 
centration (mass/volume) or as a mass flux (mass/time). 

14. Other quantities which the modeler must specify are the number of 

Stations in the longitudinal direction, Nxs, and the distance, delx, between 
them; as well as, the number of stations in the lateral direction, Nzs, and the 


lateral distance between then, delz. These parameters mav be selected as the 


modeler chooses, and in many cases the modeler may elect to make several runs 


with different selections of delx and delz. 
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15. Data for the model is setup by a program called FR’T. On typing the 
command FRNT into the microcomputer, the system responds with a series of ques- 


tions: 


STREAM VELOCITY (M/SEC) = ? 

DEPTH IN METERS = ? 

FRICTION FACTOR (RANGE = 0.005 TO 0.025) = ? 
LONGITUDINAL DISTANCE BETWEEN STATIONS IN METERS = ? 
LATERAL DIST. BETWEEN STATIONS IN METERS = ? 
NUMBER OF STATIONS IN LONGITUDINAL DIRECTION 
NUMBER OF STATIONS IN THE VERTICAL DIRECTION 
NUMBER OF STATIONS IN THE LATERAL DIRECTION = ? 
WIDTH OF INFLOW WINDOW (METERS) = ? 

STATION NUMBER FOR BOTTOM OF WINDOW (1 IS BOTTOM) 
STATION NUMBER FOR TOP OF WINDOW (< NUMBER OF Ys) = ? 
MASS FLUX OR CONC. BOUNDARY CONDITION ( M OR C) = 
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? 


‘— 


if response = M 
then: 
MASS FLUX (IN UNITS OF GRAMS/SEC) = 
else if response = C 
then: 


CONCENTRATION (GRAMS/LITER) = ? 

NUMBER OF SIZE FRACTIONS (1 to 9) = ? 

% OF MASS FOR FRACTION 1 ( O TO 100) = ? 
SETTLING VELOCITY FOR FRACTION 1 (CM/SEC) = 


repeat until data is obtained for the number of fractions specified. 


Each of the questions is answered by inputting the appropriate information 


followed by a carriage return. FRNT generates a file called PLUME.VRS (plume- 

variables). This file looks like: 
l. U = 0.05 
2. #H = 100 
3 FF = 0.01 
4 DELX = 1000 
5. DELZ = 100 
6. NXS = 41 
7 NYS = 5 
8. NZS = 5 
9. 2*BO = 50 
10. WBOT > 3 
ll. WTOP - 4 
12. OPT = M 
13. FLUX = 100000 
14. NFACS = l 
15. PCTMF = 100 
16. V - 0 

ll 


Entries in the last column are responses which the modeler gave to FRNT. The 
order of these entries is the same as that of the questions which FRNT asks. 
The file PLUME.VRS may be edited, copied to the printer, etc. In editing the 


file, the last column must be properly right justified. 


Running the Model 


16. To run the plume model, the modeler simply types in MODEL. The com- 
puter will first output the relevant input variables to the printer, and then 
respond with a series of questions which are as follows: 


PRINT TABLE OF MASS FLUXS (Y OR N)? 
PRINT SUM OF CO TABLE (Y OR N)? 
PRINT TABLE OF CD VALUES (Y OR N)? 
PRINT CONCENTRATION TABLE (Y OR N)? 


To generate a table or concentrations 
ENTER THE DEPTH AT WHICH CONCENTRATION IS WANTED 
10 stop generating concentration tables 


ENTER A NEGATIVE NUMBER 
(DEPTHS ARE ROUNDED TO THE NEAREST NODAL EVALUATION) 


PRINT TABEL OF DEPOSITION RATES (Y OR 4)? 


The tables corresponding to the second and third questions are the solutions 
(if one fraction is specified) of Equations 8 and 9 respectively. If more 
than 1 fraction is specified, the second table is the solution of Equation 8 
summed over all of the fractions. The table corresponding to the first ques- 
tion gives the total mass flux and the mass flux for the various fractions 
across planes normal to the x-axis at various downstream locations. Concen- 
trations at various depths are generated in response to the fourth query, and 
as many of these tables as are desired may be yenerated. Entry of a negative 
number halts the generation of these tables. The last table which mav be 
senerated is a table of deposition rates. This table zives deposition rates 


as mass per unit area per unit time. 


Special Considerations 


17. Referring back to Equation 13 in which vertical diffusion is ig- 
nored, observe that in the time it takes a particle to advance a distance dx , 


the distance between two stations in the longitudinal direction, it will fall 
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a distance 
dy = w/U dx (17) 


Equation 14 uses two succesive grid elevations under the assumption that the 


characteristic line 


vewx /U (18) 


falls in the interval v(j) to y(j+l) . This will happen only if the dis- 
tance dh between two vertical stations is greater than or equal to dy. The 


inequality consistent with the assumption is 


dh > dy = w/U dx (19) 


If this inequality is violated the numerical solution of Equation 8 becomes 
unstable (wild unpredictable numerical output is generated). Since the depth 
h , the settling velocity w , and the mean current velocity are fixed, and 


since 


4 


dh = h / (Nys-1) ’ (20) 


poses a restriction on the values of dx and Nys which may be used. The 


most restrictive case in a given situation will be for the material with the 
greatest settling velocity. The computer checks the inequality and when it is 
not satisfied prints the mess. _e 

Dv/Dx MUST BE > OR = w/U FOR STABILITY 

WIU © gencsecess DELT/ORLA @ caneursacs 

DECREASE DELX OR NYS ... COMPUTATION ABORTED 
In a number of instances it has been found desirable to run the model first 
with a relatively small dx until the sand fraction settles out. After re- 
moving the sand fraction from the distribution the model is rerun with a cor- 
respondingly larger dx value. 

18. Suppose that there is one fraction with zero settling velocity. 
Material must be conserved within the flow and the mass flux across anv section 
normal to the x-axis will be 

mass-flux = (2boh) avg (C)U (21) 


This equation is satisfied by the numerical model if the equation 


13 


| where Nvs is the number of stations in the vertical direction, Equation 19 im- 


Nvs - l 
avg(C_) = > C (22) 


n = 2 


is used to compute the average of C . In many instances it is the mass flux 
and not the initial concentration which is known. The numerical mcedel uses 
Equations 21 and 22 to compute the upstream concentration when the mass flux 
is given. Using the divergence theorem it is not difficult to show that the 


mass flux per unit area across a horizontal surface is 
=wC+eE ck (23) 
4 y ‘'y ies 
At the bottom the boundary condition is 
3C/3y = O (24) 


Hence, the rate of deposition is just 


q = wC (25) 


PART IV: SUMMARY AND CONCLUSIONS 


Summar 


19. The plume model developed in this study incorporates the major phys- 
ical mechanisms which influence the behavior of a suspended material. However, 
the model is still simple enough to be run on relatively small computers. The 
model does, however, have limited applicability because of the many assumptions 
made, e.g. the current is considered independent of time and is constant in 
magnitude and direction. 

20. The model has been run on both 8 bit micro and 32 bit super micro- 
computers. The modeling system consists of two programs, FRNT which sets the 
problem up, and MODEL which performs the actual computations. Both of the 
programs are interactive and use a question and response mechanism to comnuni- 
cate with the user. Variables such as the current velocity are fixed by phys- 
feal cireumctsonccs while other variables sucii ad cite umber of sCatilons in 
the x, y , and z directions are free for the modeler to choose. The onlv 
restriction on the choice of these variables is the stability criterion in- 
dicated in PART III of this report. The computer code (MODEL) checks Equa- 
tion 19 and aborts the run if the criterion is not met. 

21. A number of runs have been made in testing the model. The coarsest 
material (and the size of the computer memory) limit the longitudinal distance 
which can be modeled because of the stability criterion. When both sand and 
clay material are present, the plume length which can be modeled gives an un- 
Satisfactory picture of the fate of the finer material. What has proved to be 
useful is to remove the coarser fraction from the distribution and repeat the 
computation with an appropriately larger distance between longitudinal sta- 
tions. This practice is acceptable because the various fractions are assumed 


to behave independently in the model. 


Conclusions 


22. Much more sophisticated models of sediment diffusion processes are 
available. However, these models require a substantial body of field infor- 
mation which must be obtained for each specific problem, and also require a 


substantial investment in manpower and large computer costs for calibration 


15 


and testing. In some instances the data simply aren't available and/or the 
cost of obtaining them and applying a sophisticated model is unwarranted. The 
model described in this report gives a realistic general picture of plume be- 
havior and should be particularly useful when the use of more sophisticated 


models is impractical. 
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Appendix A: Listinz of Model 
i: INTEGER FORMFD,.LF 
2: INTEGERS1 OFT, YORN 
3: DOUPLE PRECISION D(21,21 
a: DIMENSION COh(61,11),C(21)-FLUK( 6,61). 
S$ DIMENSION VCOER(61.6) /-NUM(10) 
é: DIMENSION DPTH(21),xD1S(61),2Z01 
73 PATA FORMPFD/1I2/,CD/44180.7,C0T 
8: 
9: 11 FORMAT(19xK.F10.0) 
10; 1 FORMAT(20Ox.F10.0) 
11: - FORMAT(20x.110) 
i fas 3 FORMAT(2O9xK,AlL) 
i3: 4 FORRMAT(’ FRACT. NUP. °,I2,° % OF MASS 
14: i ‘FALL VELOCITY = °,F10.5 
iS; 
| 16; DO 10 WN = 1,10 
17: id NUM (CN) =N 
18: C 
19: WRKITE(2,6030) FORMFI 
| 20: C 
213 CALL OFEN (6,11 HFLUME VES. 
aa READ(6,11)U 
a3: REALI(6,1)4 
74: READ (6,1)FF 
| 233 WAITE(D-S)UsFR.H 
26: © FORMAT(’ U = °,F10.5.° FR = ° oF 
7% REA (6,1) 0ELx 
6: READ ( 6,1)0ELZ 
=93 WRITECO,6)0EL x DELZ 
| 3¢ 6 FORMA’ DELx oF10.3-° DELS = 
31 
32 REA (6.,2)NKS 
32 READ (6, 2)NYS 
| 24: READ (6.,2)N25 
3° WRITE CO, “INKS SNYS.NZS 
36 7 ‘FORMAT NXS = +14, NYS = 14, 
2> 
; 33 REAI (6.1) 7WORS 
39 BO = TWOFO/. 
ac BO = BC 
4) WRITE(D,8) 7WORC 
4° KEAL (6,2) 18OTw 
j 47: WAITE(2,9) I1BOTw 
44: READ (6.2) 1 7TOF w 
4°. WRITECO,I3) I 7ORW 
i a 8 FORMAT(’ WIDTH OF WINDOW = °,F 10.3 
4°; 9 FORMA’ y NODE FOTTOM OF WINDOw 
4B: 13 FORMAT < Y NODE AT TOF OF WINDOw 
49: 
| 50: READ (6,3)0F 7 
Si: WRITE(2,15)0F7T,0F7 
Seas 15 FORMAT(’ OFT = °,A1,2K.13) 
$3: IF (OFT .NE. 77) GO TO 2C 
=4: REAIT 6.1 )FLURO 
ss: FNYS=nYS-1 
56: DEL HeH/ENYS 
$7; IF CIFOTW.EQ.1)1BO%W22 
se: IF CITOPW.EGQ.NYS)ITOFWENYS-1 
39: Mw= FLOATCITOPW-IBOTW)/FNYS 
60: IF (HW.LT.O.S5 )HWel 
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a= TWOROSHW 


CONCO=FLUXO/ (USAW 


WRITE(>.19>CONCO 


9 FORMAT: CONC. 


GO ¥O 3 


a’ 
FNWYS=FLOATCNYS-2) 
RATIO = FLOAT( I TOP W-IBOT 
CONC = FLUXO/(USTWOROSHSFA 


20 READ (6.1) CONT 
30 CONTINUE 
aDIS«i)=<0. 
PO 32 I = 2ewxS 
Se EPIS¢Jo exbIsi] 
2PIS«(1)20. 
PO 34 |] = 2eNZS 
3400 Oo ZHPIS( Tp eZbIS6 I 
FNYSI= NYS-1 
FwySos wrS-. 
DELHEM/FNYS! 
DEL Y=DELH 
DF THINYS)= 
sNnY 
PO 36 1 & OeNY 
J=J-1 
De 7 . eDt 7 J+ 
3e 40=— CONTINUE 
WAITE*O.21 CONC 
-i FORMA’ GAOS s 
REAL 6.2) NFRAC 
WRITE CD. 23)NFRA 
23) FORMAT( 6NUMPES 
NO 99 Ie1,Nr¢s 
TO 99 JF ,NES 
99 CoOrcl.,yreC 
10 1000 L=1.NFRA 
READ (6.1)FCNT 
READ Sel )ViL 
VOLISVOEL 10 
CONVEFSION FRM Crs 


NYSI#NYS=-1 
DC ~ I =], 
FO? 


US BOUNDARY = 
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NTRATIC 


121: 54 CO‘. I.1)=CONCFR 

ies 60 CONTINUE 

323: C 

124: C234567 

125: VRATIO=SABS( VIL) 970 

icé: FwySi = NYS-1 

1273 DELY = H/FNYS) 

1°28: GRATIO=DELY/DELX 

icy: IF (GRATIO.GE.VRATIO)GO ™O 7¢ 
130: WRITE(2,61 )VRATIO,GRATIO 
133: 61 FORMAT(’ DELY/DELx MUST FE OR = v/U FOR ST%«ILITY 
i32: i ‘ ysU = “oF10.S6 DEL Y/DELZ = oF 10.56 
i33: > * PECREASE DELx OF wYS... COMPUTATION ABORTEL 
134: STOF 

i3s: 70 CONTINUE 

i3eé: USTR=eSORT(FA/B. 2 BU 

i137: Fxed.48uUSTF 

136: FAYSIeNYS-1 

139: FaySrenrS-- 

140: FaySte1./FeySe 

141: PEL yYensFayS! 

142: MPELY*#O.Serecy 

141: vPveO.SeviL)/DELY 

144: Unxeu/DELZ 

14°: Shri} rel yee. 

146: DPyi*e1./DELcyY 

147: RATIOS (VI(LISDEL® 
146: C 

149: fO 6O I+i.Nrs 

isC PO SO Jtlen’ 

oP DileJrec 

ise Bc CONTINUE 

123: C 

i134 wrSienys- 

1sS: C 

ise Pihedove-DY. 

is? Pi ag,2oehvi 

o - . 

159 10 100 NeoeNwYS! 

16c Ye(N-])@0DEL” 

163 YMereeDEL * 

16. yeeveHDE. Y 

163: Emefxsusree 1. <The 
164: Fe eFxauartes 1,.°-YF/* 
16% E=0.S58(EmeEF 

166 DIN Wel) e-EmeSiri 
167: PON WISeUDKeD. BESS TI 
168: Din Wel) =-EF @Siryl 
169 100 CONTINUE 


171; DONYS NYS) #VCL)-EF Sb: 


17. DINYS NYSLIe¢EP sir! 

173: C 

174: UraIe1./UD" 

173: Dtle peUnu let lel 

176: DCL DeUNKIT SHC ie. 

178: pO 200 I*2.NYSi 

179: OC LIL eu Teoil.t-i 
120: Pil, eunxienpcl.,1? 


igi: DCI, IT¢lL SUDKISH(I,I+1) 
182: -00 CONT INUE 


164: 00 300 I=2,NYS 
16s: Dil.I-iL =f 
i8é: Dil, ,lyebil, 
i167: 300 CONTINUE 
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ie 8) 

@m 
noe 


iss 

igc Ne] 

1913: Surm=0.C 

i9g2: C 

193: PbO 350 Jel .nYSi 

194: 350 SUM=SUM+CO( Jen) 

19%: FLUK(L N)=2. SUSEOsHSSURBEN YTS? 
igé: Cc 


197; PO 800 N=2,NxKS 
192: F(il)=0. 


igs: PO 400 J=l.nwySi 

00: Bic IT) eCO( IT W-1)-KATIOS(COCITO1.N-12)-COCITN-) 
203 aoc CONTINUE 

an *% © — 

203: FiNYS)=0. 

204: PO S00 Tel.enys 

0S: BCT eB Jo-0(3 3-1 98Ri 1-1 
A: sor CONTINUE 

-08: B(NYSIPSBRCNYS)/DCNYS NYS 

0° LenS 

21¢ [0 600 AK=2.NYS 

213: 

2-14: C 

eis: [00 700 Telenwrt 

-16: 70 COC I eW)ekC! 

“~« Be r 

18: VCO NeLd@e-VI(LISBLC le 

719: C VCOR MAY BE THOUGHT OF AS HAVING THE DIMENSIONS 
seW C GRAMS /SQ-mETER FEF SECOND IF OWE WILL AELGARI 
22i: C Ch AS DIMENSIONLESS 

~a A 

sees Sum=0.0 

224: DO 750 J=>.NYS! 

225: SUM=SUR+CO( JN 

wa6s =¢ CONTINUE 

“aaye Cc 

--8: FLUK(L ND] 2. @USKOSHESUMBE NYS: 


-30: B00 CONTINUE 


stes WO B20 I*i.wys 

me ® FB PO 8620 JeieNxKS 

-34: COT(T JeCOTC Le sreCOiles 

osu: B20 CONTINUE 

~36: C 

me Ps 1000 CONTINUE 

-38: C 

239: WEITE(3,681> 

740: @Ri FORMAT(’ PRINT TABLE OF MASS FLUX Y OF WN ’ 
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242: 882 


248: 

249: 6883 
250: 507 
251i: 508 
252: S501 
253: 6884 
254: 

255: 

256: 

257: 

258: 870 
259: C 
260: 

261: 687 
262: 

263: 890 
264: B95 
265: C 
266: C2345 
267: 

-68: 8691 
269: 

270: 

271: 

a7: 

275: 951 
276: 

277: 

2-78: 699 
279: C 
280: 

281: 931 


READ( 3,882) YORN 

FORMAT(A1) 

IF(YORN .EQ. 78)GO TO 895 

WRITE(2,6030)FORMFD 

WRITE(2,883) 

WRITE(2,507) 

WRITE(2,508) 

WRITE(2,501)(NUM(I),I=1,NFRACS) 

FORMAT(SX,’ TOTAL MASS FLUX AND FLUX FOR FRACTIONS *) 


FORMAT (SX, ’ GRAMS/SEC ‘/) 
FORMAT ( ’ x TOTAL FRACTION NUMBER “ ) 
FORMAT(’ METERS FLUX 15,5110) 


FORMAT(’ X(METERS) ’) 
DO 890 N=1,NXS 
SUM=O. 
DO 870 L=1,NFRACS 
SUM=SUM4FLUX(L,N) 
CONTINUE 


WRITE(2,687)XDIS(N),SUM, (FLUX(L,N)»L=1,NFRACS) 
FORMAT(1X,F8.0,1FE10.2,1PE10.2,1FE10.2,1FE10.2,1FE10.2, 
i 1PE10.2,1FE10.2,1FE10.2,1FE10.2,1FE10.2,1FE10.2) 
CONTINUE 
CONTINUE 


67 


WRITE(3,891) 

FORMAT(’ PRINT SUM OF CO TABLE -Y OR WN ae ) 

READ (3,882) YORN 

IF(CYORN .EQ. 78)G0 TO 955 

JDEL=(NYS-1)74 

JYS=14JDEL 84 

IF (JYS .EQ. NYS) GO TO 899 

WRITE(3,951) 

FORMAT(’ NO. OF Ys MUST EQUAL 5,9,13... TO PRINT THIS TABLE 
GU TO 955 


WRITE (2,6030)F ORMFD 


WRITE(2,931) 

FORMAT(20X,’ TABLE OF CO VALUES (MILLIGRAMS/LITER) ‘7) 
WRITE( 2,932) (DETH(I),I=1,NYS,JDEL) 

FORMAT(’ DEPTH = ’,F7.0,4F10.0) 


WRITE(2,884) 

DO 950 I =1,NxXS 
WRITECZ,9S3)XDIS(1),(COT( Je I) e- J=1 NYS, JDEL) 
FORMAT(1X,F8.0,1FE10.2,1FPE10.2,1FE10.2,1FE10.2,1FE10.2) 
CONTINUE 

CONTINUE 

WRITE (3,957) 

FORMAT (’ COMPUTING CI, PLEASE WAIT ‘//) 


Ch(i,1)=1.0 
DO 1050 I = 2,NZS 
Chi I,I)=0.0 


EZ=0.2*HSUSTR 
DO 1200 N=2,NxXS 
X*@(N~-1) SDELX 
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301: 
302: 
303: 
304: 
305: 
306: 
307: 
308: 
309: 
310: 
311: 
312: 
313: 
314: 
315: 
316: 
317: 
318: 
319: 
320: 
321: 
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323: 
324; 
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326: 
Se-s 


328: 
329: 
330: 
331: 
332: 
333: 
334: 
335: 
336: 
337: 
338: 
339: 
340: 
341: 
342: 
343: 
344; 
345: 
346: 
347: 
348: 
349: 
350: 
351: 


352: 


353: 
54: 


cee 
wwe 


356: 
357; 

38; 
359: 
360: 


Cc 


Cc 


441 


1100 


1250 
1300 
6030 


1301 


1302 


1303 


1304 


1340 
1380 


AL =SQRT(U/( 48XEZ)) 


DO 1100 J=#1,NZS 
Z=(J-1)SDELZ 
ARG1=AL8( BO+Z) 
ARG2=AL8( BO-Z) 
ERF i =ERF (ARG1) 
TF (ARG2 «LT. O.) GO TO 441 
ERF 2=ERF (ARGO) 
CD(N,J) #0. 58 (ERFI+ERF2) 
GO TO 1100 
ERF 2=ERF (-ARG2) 
CU(CN,J)20.S58(ERF1-ERF2) 
CONTINUE 


CONTINUE 

WRITE(3,1201) 

FORMAT(’ PRINT TABLE OF CD VALUES <Y OR N ? “) 
READ( 3,882) YORN 

IF(YORN.EQ@.78)GO TO 1300 
WRITE(2,6030)FORMFD 

WRITE(2,1202) 

FORMAT(20XK,’ TABLE OF CDT VALUES ’) 
WRITE (2,1203)(ZDIS(1I),I=1,NZS) 
FORMAT(’ Z# ’,11F10.0) 

WRITE(2,884) 

oC 1250 I = 1,NxXS 
WRITE(2,887)XDIS(I),(COCI,J),J=1,nZS) 
CONTINUE 

CONTINUE 

FORMAT (A4) 


WRITE (3,1301) 
FORMAT(’ TO GENERATE A TABLE OF CONCENTRATTION: ‘7 
‘ ENTER THE DEFTH AT WHICH CONCENTRATION IS WANTED 
‘ TO STOF GENERATING CONCENTRATION TABLES: ‘/ 
‘ ENTER A NEGATIVE NUMBER °// 
‘ <DEFPTHO Ake ROUNDED TO THE NEAREST NODAL ELEVATION 
CONTINUE 
WRITE(3,1302) 
FORMAT(’ DEFTH IN METERS = ? “) 
READ(3,1303) DEPTH 
FORMAT(F10.0) 
IFCDEFTH .LT. O0.9GO TO 1400 
WRITE (2,6030)F ORMFLI 
IY = DEPTH/DELH 
REM = DEPTH -FLOAT(IY) @DELH 
IF (CREM.GT.O.S@DELH)IY=IY41 
DEF TH = FLOAT( IY) SDELH 
WAITE(2,1304) 
FORMAT(2OX,’ TABLE OF CONCENTRATIONS (MILLIGRAMS/LITER) °) 
WRITE(2,1305) DEF TH 
WRITE(2,1203)(ZDIS(1I),I=1,.N2S) 
WRITE( 2,884) 
Iy = wyYS - ly 
DO 1380 N = 1,NxS 
DO 1340 I = 1,NZS 
C(I) = COTCIY,N)S@CDCN,I) 
WRITE(2,687)XDIS(N),(C(I),I=1,NZS) 
CONTINUE 
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367: 
368: 
369: 
370: 
371: 
372: 
373: 
374: 
375: 
376: 
377: 
378: 
379: 
380: 
381: 
382: 
383: 
384: 
385: 
386: 
387: 
388: 
389: 
390: 
391: 
392: 
393: 
394: 
395: 
396: 
397: 
398: 
399; 
400: 
401: 
402: 
403: 
404: 
405: 


Cc 
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1305 FORMAT(25X,‘ AT A DEPTH OF ’,FS.0,’ METERS ‘/) 
GO TO 1320 
1400 CONTINUE 


WRITE(3,1401) 

1401 FORMAT(’ PRINT TABLE OF DEPOSITION RATES <Y OR N- ? ”) 
READ( 3,882) YORN 
IF (YORN.EQ@.78)GO TO 1500 
WRITE(2,6030)FORMFD 
WRITE(2,1402) 

1402 FORMAT(20X,°’ TABLE OF DEPOSITION RATES “) 
WRITE(2,1403) 

1403 FORMAT( 20x, ’ (GRAMS/SEC PER M-SO )’/) 
WRITE(2,1203)(ZDIS(J),J=1,NZS) 
WRITE( 2,864) 
DO 1460 I = 1,NxXS 
SUM =O. 
DO 1420 K = 1,NFRACS 
SUM = SUM+#+VCOB(I,K) 

1420 CONTINUE 
DO 1440 J = 1,NZS 

1440 C(J)=SUMSCD(I, J) 

IN THIS CONTEXT C(4U) IS A CONVENIENT DUMMY 
WRITE(2,687)XDIS(1I),(C(J),J=1,NZS) 

1460 CONTINUE 

1500 CONTINUE 


STOF 
END 


FUNCTION ERF(X) 


DOUBLE PRECISION A1l,A2,A3,A4,A5,T.F 
DATA Al»ALLAS /.254829592,-.284496736,1.421413741/ 
DATA A4,AS 4~-1.453152027,1.061405429/ 


EXFF=0.0 
IF (xX .GT. 4.2) GO TO 10 
T=1./(1.4-32759118Xx) 
P=((((ASST+AS4) BT+AS) BT+A2) BT+A1) OT 
EXPFeEXP(-xXex) 

10 ERF=1.-PaeExrr 


RETURN 
END 
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Appendix 8: Listing of FRNT 


PROGRAM FRNT 
INTEGERS1 FLDC2O) 
INTEGERS1 CRK,LF,ZERO 
INTEGERS1 YN 

Ck=13 

LF=10 

ZERO = 48 


CALL OFPEN(7,11HFLUME VAS,1) 


WRITE(3,3) 

CALL GETFLDCFLD) 

FORMAT(’ STREAM VELOCITY (M/SEC) = 7? ’) 
WRITE(7,4)(FLD(I),1=1,20) 


WRITE(3,7) 

FORMAT(’ DEPTH IN METERS = 7? “) 
CALL GETFLICFLD)D 
WRITE(7,6)LF,(CFLOCI),1=1,20) 


WAITE(3,9) 

FORMAT( ’° FRICTION FACTOR (RANGE = 0.0005 TO 0.025) =? ’) 
CALL GETFLICFLD) 

WARITE(C7,10)LF,CFLDOCI),I#1,20) 


WRITE(3,11) 

FORMAT(’ LONGITUDIONAL DIST. BETWEEN CvaTIONC ye wevercc - 9 
CALL GETFLICFLE) 

WRITE(7,12°LF CFLOCI), 321,20) 


WRITE(3,13) 

FORMAT(’ LATERAL DIST. BETWEEN STATIONS IN METERS = 7? “*) 
CALL GETFLICFLE) 

WRITE(7,14)LF -CFLDOCI),F=1,20) 


WRITE(3,15) 

FORMAT(’ NUMBER OF STATIONS IN LONGITUDIONAL DIRECTION = 7 ’) 
CALL GETFLOCFLE) 

WRITE(7,16)LF,CFLOC I), T#1,20) 


WRITE(3,17) 

FORMAT(’ NUMBER OF STATIONS IN THE VERTICAL DIRECTION = ? °) 
CALL GETFLDCFLD) 

WRITE(7,16)LF CFL ICI), 121,20) 


WRITE(3,19) 

FORMAT(’ NUMBER OF STATIONS IN THE LATERAL DIRECTION = ? ') 
CALL GETFLOCFLD) 

WRITE(7,24)LF,CFLOCI),I2#1,20) 


WRITE(3,27) 

FORMAT(’ WIDTH OF INFLOW WINDOW (METERS) = 7? *) 
CALL GETFLDCFLD) 

WRITE(7,36)LF (FLCC I), 121,20) 


WRITE(3,29) 

FORMAT(’ STATION NUMBER FOR BOTTOM OF WINDOW ¢ 1 IS BOTTOM) = 7? *) 
CALL GETFLDCFLEI) 

WRITE(7,34)LF,(FLDCI),I#1,20) 
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WRITE(3,31) 

FORMAT(’ STATION NUMBER FOR TOF OF WINDOW ( =: NUMEBE. OF Ys) = 
CALL GETFLDCFLE) 

WRITE(7,36)LF,(FLD(I),I=1,20) 


WRITE(3,33) 

FORMAT(’ MASS FLUX ORK CONC. BOUN_CONT (M OR C) 7? ’) 
CALL GETFLICFLD) 

WRITE(7,42)LF,¢(FLD(I),I1=1,20) 

VN = FLIC20) 


IF(VN .E@. 77)GO0 TO 80 

WRITE(3,35) 

FORMAT( “’ CONCENTRATION (GRAMS/LITER) = ? “) 
CALL GETFLDC(FLD) 
WRITE(7,44)LF,(FLDC(I),I1=1,20) 


GO TO 90 
CONTINUE 


WRITE(3,21) 

FORMAT(’ MASS FLUX (IN UNITS OF GRAMS/SEC) 
CALL GETFLOCFLD) 
WRITE(7,26)LF,(FLOCI),I1=1,20) 


" 
v 
~ 


CONTINUE 


WRITE(3,23) 

FORMAT(’ NUMBER OF SEDIMENT SIZE FRACTIONS (1 70 9) = 7 °) 
CALL GETFLICFLE) 

WRITE(7,26)LF,(FLICI),I=1,20) 


NUM = FLD(20)-ZERO 


DO 100 NN = 1,NUM 

WRITE( 3,25) NN 

FORMAT(’ X OF MASS FOR FRACTION’,I2,° (0 TO 100) = 7 °) 
CALL GETFLOCFLE) 

WKITE(7,32)LF,(FLDCI),IT=1,20) 


WRITE( 3.5) NN 

FORMAT(’ SETTLING VELOCITY FORK FRACTION’,I2,’ (CM/SEC) = 7 ") 
CALL GETFLICFLD) 

WRITE(7,-6)LF,(FLECI),IT=1,20) 


CONTINUE 
12342678 

FORMAT(’ U = ’,20A1) 
FORMAT(AI, ’ Vo = ’,20A1) 
FORMAT(A1,’ # = ’,20A1) 
FORMAT(A1L,’ FF = ’,20A1) 
FORMAT(A1,’ DELX = ‘’,20A1) 
FORMAT(A1,’ DEL Z = ’,20A1) 
FORMAT(A1L,’ NXS = ‘,20A1) 
FORMAT(AL.’ NYS = ’,20A1) 
FORMAT(A1,’ NZS = ’°,20A1) 
FORMAT(A1L,’ FLUX = ’,20A1) 
FORMAT(A1,’ NFRCS = ’,20A1) 
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20 
40 


20 


= *,20A1) 
= *,20A1) 
FORMAT(AL,”’ WTOF = *,20A1) 
z 
= 
= 


FORMAT(AL,’ PCTMF 
FORMAT(A1,’ WBOT 


*,20A1) 
“»20Al1) 
*»20Al) 


FORMAT(A1L,’ 2850 
FORMAT(AL,’ OFT 
FORMAT(AL,’ CONC 
FORMAT(AL) 


WAITE(7,46)LF 
STOF 
END 
SUBROUTINE GETFLIOCFLD) 
INTEGERS] CHARS(80),.FLD(20),SF 
READ( 3.1) (CHARS(I),1#1,80) 
FORMAT(A1,B0A1) 
NCHARS = 60 
CALL COUNT(CHARS,.NCHARS) 
CALL FIELD(CHARS.FLD.NCHARS) 
RETURN 
END 
SUBROUTINE COUNT(CHARS,NCHARS) 
INTEGERS1 CHARS(80).SF 
SF = 32 
NCHARS = B80 
pO 20 I = 1,80 
NWCHARS = WCHARS - 1 
If (CHARS (NCHARS).NE.SFE GO TO 40 
CONTINUE 
CONTINUE 
KE TURN 
END 
SUPROUTINE FIELD (CRASS, FLI.NCHARS) 
INTEGERS: CHARS(60),FLI(20) 
INTEGERS1 SP 
SF = 32 
NSPS = 20 - NCHARS 
DO 20 I = 1,NSFS 
FLOCI) = SF 
CONTINUE 
PO 40 I * 1,NCHARS 
N= 1 @ WSFS 
FLICN) © CHARS(I) 
CONTINUE 
KE TURN 
END 
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